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ABSTRACT 


We present here a field geochemical study 
of controls on carbonate weathering within 
rapidly circulating, shallow groundwater— 
surface water systems in the glaciated mid- 
continent region. Groundwaters and _ sur- 
face waters in three watersheds spanning 
the Upper to Lower Peninsulas of Michigan 
consist of Ca**-Mg”-HCO, solutions derived 
from the open-system dissolution of calcite 
and dolomite in soils developed on mixed min- 
eralogy glacial drift. The thermodynamic sta- 
bilities of calcite and dolomite both decrease 
with decreasing temperature, with dolomite 
more strongly affected. Thus, the low mean 
annual temperature of these temperate 
weathering environments maximizes the 
absolute solubility of dolomite as well as its 
solubility relative to calcite. Many groundwa- 
ters in the study area approach equilibrium 
with respect to the more soluble dolomite and 
are moderately supersaturated with respect 
to calcite. Groundwaters in each watershed 
have distinct and relatively narrow ranges of 
carbon dioxide partial pressure (Poo,) values, 
which increase significantly from north to 
south (log Poo, of -3.0 to —2.2 atm), suggesting 
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that there are landscape-level differences in 
carbon transformation rates in soil weather- 
ing zones. Increases in weathering-zone FP, 
values produce HCO, concentrations that 
vary by a factor of five, but the Mg?*/Ca* 
and Mg”/HCO, ratios of all groundwaters 
are similar, suggesting relatively constant 
weathering input ratios of calcite and dolo- 
mite. Although surface waters commonly 
are between 2 and 10 times supersaturated 
with respect to calcite, the Mg”*/HCO, ratios 
of surface waters are very close to initial 
groundwater values, suggesting that back 
precipitation of calcite is not a significant pro- 
cess in these systems. The enhanced solubility 
of dolomite at low temperatures coupled with 
the landscape-level differences in carbon 
cycling suggest that temperate-zone weather- 
ing reactions in glaciated terrains are signifi- 
cant contributors to continent-scale fluxes of 
both Mg” and HCO,. 


Keywords: carbon, rivers, groundwaters, Mich- 
igan, carbonate saturation, dolomite dissolution, 
carbon dioxide. 


INTRODUCTION 


Weathering of carbonate and silicate miner- 
als in the terrestrial near-surface environment 
plays a major role in regulating atmosphere and 
ocean chemistry on time scales of millions of 
years (e.g., Holland, 1978; Berner et al., 1983). 
Much of the weathering occurs in soils at the 
interface between the biosphere and the litho- 
sphere, the site of rapid terrestrial organic matter 


production and degradation. The global riverine 
fluxes of Ca**, Mg”, and HCO, are derived 
from both carbonate and silicate sources, and 
the relative contributions from each source are 
important for the long-term CO, cycle and ocean 
chemistry. On long time scales, Mg** added to 
the oceans via carbonate and silicate weather- 
ing may be removed by exchange for Ca** via 
basalt-seawater reaction (e.g., Hardie, 1996) or 
by dolomite formation (e.g., Holland and Zim- 
mermann, 2000). Because the Mg”* supplied to 
the oceans by rivers is not removed entirely as 
a carbonate, the input fluxes of HCO, from the 
weathering of Mg-carbonates and Mg-silicates 
are not balanced by equivalent output fluxes of 
Mg-carbonates. Thus, the relative riverine inputs 
of Mg” from carbonate versus silicate sources 
bear significantly on the global carbon cycle. 
On shorter time scales, the predicted response 
of terrestrial systems to global warming and 
atmospheric CO, increases is expected to include 
increases in biomass and in soil respiration rates. 
The terrestrial carbon reservoir has commanded 
increasing attention as the likely sink for anthro- 
pogenic CO, (Houghton et al., 1998, 1999). An 
increase in terrestrial biomass may result from 
the “fertilization” effect of elevated CO, levels 
on photosynthesis, with a large portion of this 
partitioned into below-ground biomass (e.g., Van 
Veen et al., 1991; Bowden et al., 1993). Plants 
grown in elevated CO, atmospheres have more 
rapid root and microbial respiration rates, which 
can increase CO, partial pressure (P.,.) values 
in soil gases (Hungate et al., 1997; Schlesinger 
et al., 2000; Andrews and Schlesinger, 2001; 
Williams et al., 2003). The increase in soil gas 
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CO, may simply lead to an enhanced flux back 
to the atmospheric reservoir (e.g., Trumbore et 
al., 1996), or it may drive an increased carbon 
flux to the soil carbon pool and/or to groundwa- 
ter systems as dissolved inorganic and organic 
carbon (e.g., Richter et al., 1995; Schlesinger, 
1997; Williams et al., 2003). 

Rapid equilibration of soil waters with car- 
bonate minerals has been observed in experi- 
mental tree-growth chambers (Williams et al., 
2003) and in natural shallow groundwaters 
(e.g., Freeze and Cherry, 1979; Reardon et al., 
1979, 1980; Seifert, 1982; Van Breenan and 
Protz, 1988; Kempe et al., 1991; Ku, 2001). It 
is thus very important to understand the links 
between the rates of mineral weathering and 
organic matter dynamics in terrestrial settings. 
Via interaction with mineral solubilization, soil 
reactions police the rate of atmospheric carbon 
transfer from mid-latitude ecosystem reservoirs 
to hydrospheric reservoirs (groundwaters and 
oceans). Increased fluxes of dissolved inorganic 
carbon (DIC) and cations from mid-latitude 
areas also are important in terms of assessing 
changes to surface water compositions and for 
inputs to the global ocean reservoir. The annual 
global flux of dissolved carbon as DIC and dis- 
solved organic carbon (DOC) from rivers has 
been estimated at ~0.6 Gt C/yr (e.g., Kempe, 
1985; Kempe et al., 1991; Aumont et al., 2001), 
which is a small, but nontrivial, contribution to 
the global carbon budget. 

Although there have been many studies of 
mineral weathering in surface environments, 
few have focused on the controls on carbon- 
ate equilibria in such systems or on the relative 
contributions of dolomite and calcite to weath- 
ering fluxes. Studies on solute fluxes from the 
terrestrial environment to the oceans via riverine 
discharge have focused primarily on the down- 
stream reaches of major world rivers such as the 
Amazon, Yellow, and Mississippi Rivers (e.g., 
summary in Meybeck, 1987). The large drainage 
basins of these rivers integrate varying climates, 
lithology, soil compositions and depths, and 
land-use types, including agriculture, urban, and 
forest. Rates of chemical weathering and their 
implications for geochemical cycling over long 
time scales have been investigated (e.g., Blum et 
al., 1998; Huh et al., 1998; Jacobson et al., 2002; 
Millot et al., 2003). Others have focused on the 
internal riverine processing of trace elements 
and nutrients, isotopes (O, C, Sr) (Palmer and 
Edmond, 1992; Yang et al., 1996; Karim and 
Veizer, 2000), or the organic/inorganic specia- 
tion of carbon (Aucour et al., 1999; Roy et al., 
1999; Telmer and Veizer, 1999; Finlay, 2003; 
Hélie et al., 2002; Barth et al., 2003). 

The current fluxes of inorganic carbon from 
each continent are available from compilations 
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of annual river-water discharge and DIC con- 
tents (e.g., Ming-hui et al., 1982; Kempe, 1985; 
Berner and Berner, 1996; Huh et al., 1998). 
Nearly 85% of the 0.6 Gt C/yr flux is derived 
from the Northern Hemisphere continents of 
Europe, North America, and Asia, all of which 
have significant areas of sedimentary bedrock 
and unconsolidated, glacially derived material 
(Holland, 1978; Meybeck, 1987). In general, 
the solubility of feldspars is very low, and dis- 
solution rates are thousands of times slower than 
those of carbonate minerals (e.g., Holland, 1978; 
Amiotte Suchet and Probst, 1993; Berner and 
Berner, 1996; Roy et al., 1999). Rivers exclu- 
sively draining continental silicate regions tend 
to have very low HCO, concentrations (i.e., 
<0.5 meq L" as HCO, ; Holland, 1978; Berner 
and Berner, 1996), whereas rivers draining car- 
bonate-bearing areas, such as the Danube River 
in Europe and the Yellow River in China, tend to 
have higher overall HCO, concentrations and 
HCO, fluxes. 

In North America, almost 90% of the HCO, 
flux is derived from the mid-continent region, 
where the bedrock geology is dominantly sedi- 
mentary and mantled by carbonate-bearing 
glacial drift. It is this mid-continent region that 
we selected for an integrated study of surface 
water and groundwater geochemistry to better 
constrain landscape-level controls on carbonate 
weathering. Figure 1 is a schematic of the Mis- 
sissippi and St. Lawrence River basin drainage, 
large portions of which are established atop land- 
scapes that have been affected by repeated cycles 
of Pleistocene glaciation. Given the rapid reac- 
tion rates of carbonate minerals, the chemistry of 
such drainage systems is likely to be quite sensi- 
tive to short-term changes in the carbon cycle. 
According to Raymond and Cole (2003), there 
appeared to be a significant increase in the alka- 
linity flux from the Mississippi drainage system 
during the past half-century, an increase they 
attribute to enhanced fluxes from agriculturally 
dominated watersheds. This is consistent with 
data on soil CO, contents and fluxes (Reardon 
et al., 1979; Brook et al., 1983; Buyanovsky and 
Wagner, 1983; Solomon and Cerling, 1987; Lee, 
1997; Andrews and Schlesinger, 2001), which 
show that soil gas CO, profiles are generally 
controlled by depth variation in rates of micro- 
bial and root respiration and diffusive transport. 

Early studies of the thermodynamic controls 
on the carbonate geochemistry of near-surface 
waters primarily focused on the role of CO, 
partial pressure on calcite solubility, rather 
than temperature, as the most important vari- 
able (e.g., Brook et al., 1975; Harmon et al., 
1975; Drake and Wigley, 1975; Trainer and 
Heath, 1976; Drake, 1983). However, another 
important control on carbon fluxes is the inverse 


temperature dependence of carbonate mineral 
thermodynamic stability. The stability of both 
dolomite and calcite decreases with decreasing 
temperature, but the temperature effect is stron- 
ger for dolomite, which has an equilibrium sta- 
bility constant that is more strongly dependent 
on temperature (e.g., Langmuir, 1997; Sher- 
man and Barak, 2000; Hyeong and Capuano, 
2001). Using the most recently compiled ther- 
modynamic data for dolomite (e.g., Hyeong and 
Capuano, 2001), the equilibrium constant for 
dolomite becomes progressively greater than 
that of calcite as temperatures decrease below 
~30°C (Fig. 2A). Considering a solution at 
10 °C and a Mg*/Ca** molar ratio of 0.5, the 
solubility of calcite versus dolomite at various 
Pee: values, in terms of HCO, total concentra- 
tion, is significantly enhanced relative to pure 
calcite (Fig. 2B; after Williams et al., 2003). 

The Michigan landscape provides an excel- 
lent natural laboratory in which to examine 
controls on groundwater and surface water 
geochemistry and the limitations they impose 
on carbonate weathering intensities. In order 
to better understand the fundamental controls 
on these fluxes, we investigated the carbonate 
mineral equilibria and dissolution-precipita- 
tion mass balances in watersheds developed on 
Pleistocene glacial drift deposits, many of which 
are carbonate-bearing. We integrated data on the 
carbonate geochemistry of surface waters and 
groundwaters in three Michigan watersheds, the 
Tahquamenon, the Cheboygan, and the Huron, 
located along an environmental gradient (loca- 
tions shown in Fig. 3). Catchment topography 
is dominated by upland moraine systems, which 
control drainage divides on both regional and 
local watershed scales. Close chemical linkages 
between groundwater and surface water sys- 
tems have been forged by the high permeability 
of soils and underlying glacial drift across the 
state (Bent, 1971; Mandle and Westjohn, 1989; 
Rheaume, 1991; Holtschlag, 1996; Grannemann 
et al., 2000). The vast majority of groundwaters 
discharge into the Great Lakes catchment, with 
<2% of meteoric water discharging into the 
deeper basinal-scale flow system (Eberts and 
George, 2000). 


HYDROGEOLOGIC FRAMEWORK OF 
FIELD STUDY SITES 


In designing our field study, we considered 
differences in glacial history, drift and bedrock 
lithology, as well as temperature, geochemi- 
cal and ecological gradients. Accordingly, we 
selected three watersheds in the state of Michigan 
for comparative investigation: the Tahquamenon 
watershed in the Upper Peninsula, the Cheboy- 
gan watershed in the northern Lower Peninsula, 


Geological Society of America Bulletin, May/June 2007 


Carbonate geochemistry in temperate glaciated watersheds 


Mississippi River and St. Lawrence River Drainage 


Figure 1. Hydrological regions of central North America showing the St. Lawrence and Mississippi River drainages. The significant regional 
extent of Pleistocene glaciation on these two major drainages is shown schematically. 
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Dolomite vs. Calcite Solubility (Mg/Ca = 0.5, T = 10 °C) 
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Figure 2. Relevant stability data for calcite versus dolomite. (A) Dolomite and calcite thermodynamic data (K,,,,,. and 2 K ) as function 
of temperature showing that the equilibrium constant value for dolomite increases more significantly with decreasing temperature than does 
the equilibrium constant value for calcite (after Langmuir, 1997). (B) Solubility at 10 °C, in terms of HCO, total concentration of pure calcite 
versus dolomite in a solution with Mg”*/Ca* molar ratio of 0.5 at various values of CO, partial pressure (after Williams et al., 2003). 


dolomite 
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A Topography 


Tahquamenon 


B Bedrock lithology 


Precambrian igneous, 
sedimentary and 
metamorphic rocks 


Sandstone 
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Figure 3. Location maps of watershed study sites in the state of Michigan with the study 
watersheds indicated by shading: the Tahquamenon (Upper Peninsula), the Cheboygan, and 
the Huron watersheds from north to south. (A) Shaded relief map, with topography gener- 
ally controlled by variation in thickness of glacial deposits. (B) Subcrop bedrock lithology 
under Pleistocene glacial deposits. “(Red Beds” includes multiple Jurassic terrestrial and 
marginal marine sedimentary deposits that consist mainly of sandstone, shale, and clay, 
with minor beds of limestone and gypsum. 


and the Huron watershed in southeastern Lower 
Peninsula. The general framework of these mid- 
continent surficial aquifer systems and drainage 
basins is shown in Figures 3A, 4A, 4B, and 4C. 
The Tahquamenon drains into Lake Superior, 
and the Cheboygan and the Huron drain into 
Lake Huron and Lake St. Claire, respectively, 
and all are part of the St. Lawrence—Great Lake 
drainage system. The significant latitudinal 
gradient represented by our study sites leads to 
differences in mean annual temperature, which 
ranges from a low of 4.8 °C in the Tahquamenon 
watershed, to 10.0 °C in the Huron watershed 
near Ann Arbor. 

Glacial drift deposits commonly contain 
carbonate minerals reworked from underlying 
bedrock formations during the multiple glacial 
advances and retreats. Soil development over 
the past 6000—12,000 yr has led to accumulation 
of organic carbon at greatly accelerated rates 
relative to nonglaciated landscapes (Harden et 
al., 1992). Largely unaltered parent material is 
typically encountered within 1-3 m of the sur- 
face. Dissolution of carbonate minerals may 
be enhanced in glacial drift deposits due to the 
grinding and stress induced by glacial actions 
(e.g., Sharp et al., 1995, 1999; Anderson et al., 
2000). Carbonate dissolution profoundly affects 
soil solution pH, soil profile development, metal 
hydroxide distributions, and solute fluxes to 
groundwaters (Reardon et al., 1980; MacDon- 
ald et al., 1991, 1992). These geochemical rela- 
tions contrast sharply with those of the more 
intensively studied carbonate-poor landscapes 
that dominate the eastern and southern regions 
of the United States (e.g., Likens and Bormann, 
1995; Richter et al., 1995; Lee, 1997). 

Like much of the adjacent mid-continent 
region, the study watersheds were developed 
atop variable thicknesses of glacial drift depos- 
ited in the waning stages of the last glacial 
advance. Thicknesses of up to 300m are not 
uncommon, and many groundwater supply 
wells in Michigan tap glacial drift aquifers (e.g., 
Rheaume, 1991). The timing of glacial retreat 
varies across the watersheds: the Huron water- 
shed was established on drift deposits ca. 18 ka 
in age, while the Cheboygan and Tahquamenon 
watersheds in the northern Lower and Upper 
Peninsula of Michigan are only ca. 11 ka in age 
(e.g., Futyma, 1981). The drift deposits over- 
lie Paleozoic bedrock units of variable lithol- 
ogy, permeability, and porosity (see Fig. 3B). 
In brief, carbonate rocks are most abundant in 
the Silurian-Devonian sections, which have sub- 
crops located along the northern and southern 
rims of the Michigan Basin. These carbonates 
are commonly interbedded with evaporite min- 
erals and consist of both calcite and dolomite 
interbeds (e.g., Catacosinos et al., 1990). Thus, 
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glacially reworked sediments throughout the 
Lower Peninsula of Michigan commonly con- 
tain both carbonate minerals, as has been dem- 
onstrated by studies of drift soil mineralogic 
profiles (Schaetzl, 1992; Jin, 2007). 

The Tahquamenon watershed has very sub- 
dued topography and largely represents slow 
drainage from wetlands established atop thin, 
quartz-rich glacial drift. The glacial drift man- 
tles sandstones of Proterozoic and Cambrian 
age and Ordovician shales that represent the 
base of the Michigan Basin stratigraphic sec- 
tion. The Cheboygan and Huron watersheds 
have much more dramatic topography due to 
thick moraine deposits in these areas. Carbonate 
mineral content of the glacial drift in the Che- 
boygan and Huron watersheds is much higher 
than in Tahquamenon, and it affects the water- 
shed solution chemistry and soil profile devel- 
opment. Results from previous hydrogeochemi- 
cal studies of Michigan surface water chemistry 
suggest that there are fundamental differences 
in surface water chemistry among these water- 
sheds (Rheaume, 1991). 

For the Tahquamenon watershed, the mean 
annual temperature is 4.8 °C, and the average 
annual precipitation is 82cm (NOAA, 1992). 
The Quaternary sediments in the Tahquame- 
non watershed are thin (<15 m) and are mainly 
composed of lacustrine organic and clay depos- 
its with some ice-contact moraines, including 
sands and silts. The underlying bedrock geol- 
ogy is dominantly shale and sandstone (see 
Fig. 4B). The watershed has an ~100 d growing 
season, and vegetation is dominated by white- 
red pine and Jack pine forests, some northern 
hardwoods, and pasture land and limited crops 
(Albert, 1995). 

The Cheboygan watershed has a mean annual 
temperature of 5.5 °C, and the average annual 
precipitation is 84cm (NOAA, 1992). Glacial 
deposits range between 15 and 180 m thick in 
the watershed. Thick deposits of highly per- 
meable lake beach sands form the floor of the 
watershed as a result of highstands of the Great 
Lakes (Farrand, 1982), while sandy glacial till 
forms the highlands. The Cheboygan watershed 
has the thickest glacial deposits of the three 
watersheds, and as a result, the greatest relief. 
The bedrock geology underlying the glacial 
deposits is Devonian-Silurian limestones and 
dolostones (see Fig. 3B). The watershed has a 
growing season of 100 d, and the dominant land 
use is northern hardwood forests and pasture 
(Albert, 1995). 

The Huron watershed has a mean annual tem- 
perature of 10.0 °C, and the average annual pre- 
cipitation is 85 cm (NOAA, 1992). The glacial 
deposits are primarily contact features including 
ground moraines and kames with some outwash. 
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Figure 4. Watershed drainage map sketches for each of the three study watersheds with 
approximate groundwater locations and surface water locations indicated. Groundwater loca- 
tions are shown in closed circles, and surface water locations are shown in open triangles. 
(A) The Tahquamenon watershed. (B) The Cheboygan watershed. (C) The Huron watershed. 


The moraines are composed of silts, sands, and 
gravels, while the outwash is mostly sorted sands 
and gravels (Farrand, 1982). The Quaternary 
sediment is variable in thickness (<15—122 m) 
and is inhomogeneous, containing lenses of 
material of differing permeability. The bedrock 


geology underling the Huron watershed is pri- 
marily Devonian shales, with sandstone along 
the eastern edge (see Fig. 4C). Studies on car- 
bonate precipitation (Szramek and Walter, 2004) 
and arsenic mobility (Szramek et al., 2004) have 
been conducted on the Huron watershed. The 
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Huron watershed has a growing season of 140- 
150 d, and oak-hickory and beech-sugar maple 
forests dominate the vegetation. The longer 
growing season allows for extensive croplands 
rather than pastures (Albert, 1995). 


METHODS 
Field Sampling 


Sampling in the three watersheds was con- 
ducted between 1997 and 2004. Groundwater 
was collected in summer 2000 and late fall 
2003 in the Tahquamenon, in fall 1997, 2000, 
and spring 2003 in the Cheboygan, and in fall 
2000, 2001, and spring 2001 and 2002 in the 
Huron watershed. Surface water samples from 
the three watersheds were collected seasonally 
(early spring through fall) to capture annual 
variations in chemistry. The Tahquamenon 
streams were sampled during 2000-2003, 
the Cheboygan in 1997-2003, and the Huron 
in 1999-2001 and in 2003. Stream sampling 
locations for each watershed are shown in 
Figures 4A-4C. At sampling sites, stream 
waters were collected upstream of large roads 
and developments to limit potential contami- 
nation from local runoff. About 240 surface 
water samples were taken in total. Groundwa- 
ter was collected from household and small 
business wells that tap shallow Quaternary 
aquifers. Data from drillers’ logs (depth, sub- 
surface geology, static water level) were used 
to guide selection of wells for sampling. About 
75 groundwater samples were obtained from 
the three watersheds. 

Temperature, conductivity, and pH for each 
sample were determined in the field. Tem- 
perature was measured using a YSI model 58 
meter and a YSI 5239 probe with high-sensi- 
tivity membranes directly at the source, either 
in the stream or at the groundwater well in the 
flow stream. Conductivity was measured using 
a Corning 316 meter to provide a rapid geo- 
chemical reference point, with a precision of 
+5%. A Corning 315 high-sensitivity pH meter 
with an Orion Ross combination pH electrode 
was calibrated with low ionic strength buffers 
of 4.10 and 6.97 (corrected for temperature) 
and used to measure pH in the field as close 
to the water temperature as possible. Because 
the pH of a sample can change due to degas- 
sing and warming, samples were placed in a 
large-volume airtight container and measured 
at least twice to ascertain electrode stability. 
The reproducibility of field pH determinations 
was +0.02 pH units. 

Sample aliquots for later chemical analysis 
in the laboratory were field-filtered through a 
0.45 um nylon filter into their respective bottles 
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and kept cold until analyzed. Acid-washed 
high-density polyethylene (HDPE) bottles 
were used to collect cation aliquots. HDPE 
bottles cleaned with purified water were used 
to collect alkalinity and anion samples. Double 
deionized water—-washed glass serum bottles 
were used to store samples for DIC analysis. 
DIC samples were poisoned in the field against 
biological activity using CuCl, and were filled 
to the top with no headspace before capping. 
Cation aliquots were preserved in the field with 
ultrapure HNO,. 


Laboratory Analyses 


Major-element chemistry was measured by 
inductively coupled plasma—atomic emission 
spectroscopy (ICP-AES) for cations and ion 
chromatography (IC) for anions on several dif- 
ferent instruments, due to laboratory equipment 
upgrades. Prior to 2002, a Leeman Labs, Inc., 
Plasma-Spec ICP-AES 2.5 was used to analyze 
Ca?*, Me", Na‘, K*, Sr**, Ba’, Mn”, Fe?**, 
Zn**,H ,Si0,, and B with a precision of +2% for 
major and +5% for minor elements. Samples 
from 2002 were analyzed on a Perkin Elmer 
Optima 3300 DV, and samples from 2003 were 
analyzed on a Jobin Yvon Horiba Ultima 2C. 
Analytical precision was +4% on the Leeman 
ICP and +2% on the Perkin Elmer and Jobin 
Yvon Horiba instruments. Anions (F, CI, 
SO,*, and NO,, F,, and Br) were analyzed 
on a Dionex 4000I series ion chromatograph 
(IC) with an AS14 column prior to 2002 and 
a Dionex ICS-2500 IC in 2002-2003. Analyti- 
cal precision was +3% on the 4000I series IC 
and +2% on the 2500 series IC. Reanalysis of 
several ICP and IC samples on the new instru- 
ments yielded results consistent with previous 
measurements, within analytical precision. Br 
and F were consistently below detection and 
were excluded from the results table. Samples 
gathered for DIC (= H,CO,° + HCO, + CO,”) 
were analyzed on a UIC Coulometrics CO, 
coulometer with a precision of +2%. Titration 
alkalinity (here largely equal to carbonate alka- 
linity, HCO, + 2CO,”) was measured by elec- 
trometric endpoint titration using a Radiometer 
TitraLab automated titration system. Due to the 
given measurement precision (+0.01 meq/kg) 
and the pH range of the samples, the carbonate 
alkalinity was equivalent to the concentration 
of HCO,. Charge balance calculations were 
always within +5%. The measured DIC and 
titration alkalinity measurements of the sam- 
ples agreed well within the precision of the two 
measurements, suggesting that the CO, system 
was well constrained and consistent among the 
three carbonate system parameters measured 
(e.g., pH, DIC, and titration alkalinity). 


Thermodynamic Modeling and Aqueous 
Speciation 


Aqueous speciation and mineral ion activ- 
ity products (IAPs) were determined using 
the U.S. Geological Survey (USGS) program 
SolminEQ.88 (Kharaka et al., 1988). The car- 
bonate system parameters used were usually car- 
bonate alkalinity and pH, although results with 
DIC and pH were cross-compared to determine 
internal consistency. Values for Pog, IAP/K it. 
(Qceal), and IAP... Were obtained from the 
SolminEQ.88 output. As discussed by Drever 
(1997), the normalized saturation state ratio con- 
cept from Zhang and Nancollas (1990) was used 
to compare the IAP of dolomite [CaMg(CO,),] 
on the same molar basis as that of calcite 
(CaCO,). This was accomplished by taking the 
square root of IAP... rom SolminEQ.88. 
This normalized IAP... Value was com- 
pared to the “log K,,,.,,. (equilibrium constant) 
value (—8.45) at 25 °C (Hyeong and Capuano, 
2001) to obtain Q,.. nie for sedimentary dolo- 
mite of intermediate ordering. Given a preci- 
sion of +0.02 pH units and the uncertainty in the 
titration alkalinity, the uncertainty in P,,/ was 
calculated to be ~10%. The error in Q for calcite 
and dolomite was within 25% as determined by 
the precisions of the carbonate input parameters 
(DIC, alkalinity, pH). 


RESULTS AND INTERPRETATIONS 


Complete geochemical databases for the 
field, laboratory, and calculated parameters for 
groundwaters and surface stream waters are 
presented in Appendixes | and 2, respectively.' 
Averaged data are presented in Table 1 for 
groundwater samples and in Table 2 for surface 
stream water samples. 


Elemental Stoichiometries and Mineral- 
Water Reactions 


The groundwaters within these shallow glacial 
drift aquifers acquire solutes from the dissolution 
of minerals in the soil zone. As shown in Fig- 
ure 5A, groundwater samples from each of the 
three watersheds are dominantly HCO,, Ca**, 
and Mg” waters, and they fall close to a 2:1 line, 
which is consistent with the ratio obtained by 
dissolving carbonate minerals with carbonic acid 
(Ca** + Mg”: 2 HCO, ). Groundwaters from all 
three watersheds have an average (Ca** + Mg**): 


'GSA Data Repository item 2007059, detailed 
groundwater and surface water geochemistry data 
from the Tahquamenon, the Cheboygan, and the 
Huron watersheds, is available on the Web at http:// 
www. geosociety.org/pubs/ft2007.htm. Requests may 
also be sent to editing @ geosociety.org. 
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TABLE 1. GROUNDWATER GEOCHEMISTRY 


Watershed Value pH 
Tahquamenon Average 8.07 
(n= 21) Median 8.08 
STD 0.32 
Cheboygan Average a a 
(n = 28) Median 7.74 
STD 0.19 
Huron Average 7.59 
(n= 22) Median 7.61 
STD 0.12 


T Alk = Ca**+Mg* Mg*/HCO, Mg*/Ca** Log PCO, (IAP/K),,, (IAP/K).,, 
(°C) = (medq/L) (mM) (atm.) 
8.7 2.42 1.28 0.18 0.53 2.96 0.96 1.42 
8.4 2.47 1.28 0.19 0.51 -2.95 0.98 1.40 
2.0 0.60 0.44 0.03 0.12 0.36 0.43 0.62 
9.0 3.54 1.90 0.21 0.64 -2.45 0.65 1.07 
8.7 3.58 1.80 0.20 0.63 —2.52 0.66 1.05 
2.5 0.90 0.50 0.03 0.08 0.30 0.17 0.29 
11.3 5.31 2.77 0.17 0.51 -2.13 1.46 1.91 
W141 5.37 2.71 0.17 0.50 -1.95 1.35 1.85 
0.9 0.30 0.31 0.01 0.05 0.13 0.39 0.45 


Note: Alk—alkalinity, |AP—ion activity product, STD—standard deviation, K—equilibrium constant, dol—dolomite, cal—calcite. 


TABLE 2. SURFACE WATER GEOCHEMISTRY 


Watershed Value 
Tahquamenon Average 
(n= 42) Median 
STD 
Cheboygan Average 
(n= 138) Median 
STD 
Huron Average 
(n= 66) Median 
STD 


pH T (°C) Alk Ca**+Mg* Mg**/HCO,- 
(meq/L) (mM) 
7.54 12.77 1.18 0.71 0.20 
7.62 11.10 1.22 0.76 0.20 
0.51 8.19 0.56 0.29 0.05 
8.17 12.36 2.84 1.54 0.20 
8.21 12.85 2.96 1.60 0.20 
0.30 6.90 0.45 0.27 0.02 
8.15 16.77 4.07 2.68 0.23 
8.20 19.65 3.96 2.62 0.22 
0.28 7.19 0.56 0.42 0.05 


Mg**/Ca** Log PCO, (IAP/K),., 
(atm.) 
0.47 —2.76 0.54 
0.46 -2.73 0.12 
0.08 0.37 0.83 
0.59 -3.00 3.27 
0.58 -3.01 2.81 
0.11 0.26 2.50 
0.52 2.22 7.37 
0.51 2.74 6.00 
0.10 1.24 5.00 


Note: Alk—alkalinity, |AP—ion activity product, STD—standard deviation, K—equilibrium constant, cal—calcite. 
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Figure 5. Ca* + Mg** versus HCO, relationships for groundwater (A) and surface streams (B) within the three study catchments. The 
groundwaters for the three watersheds fall close to the 2:1 line (2HCO, = Ca + Mg”), and no significant difference between the ground- 
waters is observed according to a (Ca** + Mg”*)/HCO, ratio of 0.51 + 0.01. This strongly indicates that the groundwater chemistry is con- 
trolled by carbonic acid dissolution of calcite and dolomite present in the glacial drift. 


HCO, ratio of 0.53 + 0.06, indicating no signifi- 
cant difference in their relative solute distribu- 
tions. Groundwaters from the Huron watershed 
have the highest alkalinity values (mean HCO, 
concentration of 5.31 +0.30 mM), whereas 
groundwaters from the Tahquamenon watershed 
have much lower alkalinity values (mean HCO, 
concentration of 2.42 + 0.60 mM), less than half 
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those of the Huron watershed. These watershed- 
scale differences suggest that carbonate mineral 
solubility within each watershed may be con- 
trolled locally by factors such as climate, abun- 
dance of carbonate, soil thickness, and physical 
aquifer characteristics. 

The major solute composition of surface 
waters in these watersheds, shown in Figure 5B, 


is influenced by contributions from groundwater, 
precipitation, runoff from anthropogenic sources, 
and in-stream or lake processes. Shallow glacial 
drift groundwater influx, however, is the major 
source of water and solutes to these surface sys- 
tems, as can be seen from the similarity between 
groundwater and surface water major-element 
chemistries (compare groundwaters and streams 
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B Streams 


Cl (mM) 
Figure 6. SO,? versus Cl anion relationships for groundwaters (gw) (A) and surface streams (B) within the three study catchments. The 


Tahquamenon watershed has higher SO,” to Cl ratios than do the Huron and Cheboygan watersheds. The Huron watershed has particu- 
larly elevated Cl concentrations in surface waters due to anthropogenic inputs from road de-icing salts. 


in Fig. 5). The relationship between HCO, and 
Ca** + Mg concentrations in surface waters 
again indicates a stoichiometry consistent with 
carbonate mineral dissolution by carbonic acid. 
The composition of the surface waters of the 
Tahquamenon and Cheboygan watersheds falls 
close to the 2:1 line, similar to their respective 
groundwaters, but surface waters are generally 
more dilute. The Huron surface waters, unlike 
the Huron groundwaters, however, lie below the 
2:1 line, indicating a relative increase in Ca** 
and/or Mg” relative to HCO,. 

This imbalance in the divalent cation to bicar- 
bonate chemistry of the surface waters in the 
Huron watershed is largely due to the input of 
anthropogenic salts, such as CaCl, and NaCl, 
which contribute Cl and possible sulfide oxida- 
tion to produce SO,*. Although groundwaters 
from each watershed contain significant amounts 
of Cl and SO,” (see Fig. 6A), these are generally 
well below 0.5 mM. The Tahquamenon ground- 
waters have higher SO,” than Cl (SO,7/CI of 
2.5 to 8), whereas the Huron and Cheboygan 
groundwaters in general have a SO,”/CI ratio of 
less than 1. The SO,* and Cl concentrations of 
the surface stream waters are also shown in Fig- 
ure 6B. For the Tahquamenon and Cheboygan 
watersheds, surface waters have a very similar 
range of SO,> and Cl concentrations to those 
of their source groundwater values. In contrast, 
the surface waters of the Huron watershed have 
markedly higher Cl and SO,> concentrations 
than the source groundwaters, with surface water 
concentrations of Cl ranging between 0.3 and 
3.0 mM and SO,” generally between 0.02 and 
1.0 mM. The majority of Huron surface waters 
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have low SO ACE ratios, much less than 1. 
The Huron watershed experiences considerably 
more anthropogenic pollution than do the other 
two watersheds, and the large input of Cl to 
surface waters is derived from de-icing salts and 
water softeners (e.g., Jones and Sroka, 1997). As 
discussed by Szramek and Walter (2004), the 
Cl from anthropogenic salts is charge balanced 
by both Na* and Ca** in the surface waters of the 
Huron watershed, which produces the excess 
divalent cation concentrations relative to HCO, 
shown in Figure 5B. 

The relative proportions of Mg and Ca?* 
in the groundwater and surface waters of the 
three watersheds are shown in Figures 7A and 
7B, respectively. While Ca” is a product of the 
weathering of both calcite and dolomite, Mg”* is 
derived mainly from dolomite in the lower Paleo- 
zoic carbonate material found in the Michigan 
Basin. The Mg?*/Ca** ratios of groundwaters 
and surface waters are similar in the three study 
watersheds. Groundwaters of the Tahquamenon, 
Cheboygan, and Huron watersheds have aver- 
age Mg”*/Ca** ratios of 0.53 + 0.12, 0.64 + 0.08, 
and 0.51 + 0.05, respectively (see summary of 
averages in Table 1). Similarly, surface waters 
had average ratios of 0.47 + 0.08, 0.59 + 0.11, 
and 0.52 + 0.10, respectively (see summary of 
averages in Table 2). The Mg”*/Ca** mole ratios 
for groundwaters and surface waters shown in 
Figure 7 always exceed a value of 0.3 and are 
generally less than a value of 0.75. 

When considering calcite and dolomite molar 
formulas normalized to CO,”, a solution Mg”*/ 
Ca** ratio of 0.33 would be attained by dissolv- 
ing nearly equal amounts by weight of calcite 


@ Tahquamenon 
O Cheboygan 
@ Huron 


Cl (mM) 


and dolomite (one mol Ca” from calcite plus 0.5 
mol Ca** from dolomite, and 0.5 mol Mg** from 
dolomite). Similarly high Mg?*/Ca”* ratios have 
been observed in surface waters and groundwa- 
ters from other carbonate weathering environ- 
ments, most recently in a very thorough study 
of solute sources and ion exchange processes of 
the Bow River in southern Alberta (Grasby and 
Hutcheon, 2000). Here, elevated Mg”*/Ca** ratios 
were interpreted as evidence for smectite ion 
exchange reactions, even though dolomite disso- 
lution is the dominant source of Mg (Grasby et 
al., 1999). In contrast, kaolinite is the dominant 
secondary aluminosilicate in soil profiles of the 
three Michigan watersheds (Williams, 2005; Wil- 
liams et al., 2007), so we dismiss ion exchange 
mechanisms involving smectite as a significant 
control on the Mg**/Ca** ratios of groundwaters 
and surface waters here. Instead, the slight dif- 
ferences observed between the watersheds could 
be related to the relative amounts of Silurian 
and Devonian carbonate rock, which have dif- 
ferent proportions of dolostone relative to lime- 
stone lithologies, contributing to the glacial drift 
deposits. In addition, the sequence in which the 
soil waters contact calcite versus dolomite in 
systems open to CO, also influences the relative 
amounts of calcite and dolomite that can be dis- 
solved at saturation (Palmer and Cherry, 1984). 


Carbonate Speciation and Carbonate 
Mineral Equilibria 


The soil zone in the shallow, porous glacial 
drift aquifers of the Tahquamenon, Cheboy- 
gan, and Huron watersheds provides an ample 
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Figure 7. Mg” versus Ca* concentrations for groundwaters (A) and surface streams (B) within the three study catchments. The dissolution 
of both dolomite and calcite is evident from the high Mg”* and Ca* concentrations and the high Mg”*/Ca** ratios for the groundwaters and 
surface stream waters. The ratio 0.55 + 0.13 indicates that about one mol of dolomite [as CaMg(CO,),] dissolves for every 1 mol of calcite 
(CaCO,) within the weathering zone of the watersheds. 


supply of CO, for weathering in an open-system 
environment. The partial pressure of CO, within 
the groundwaters is at or below atmospheric 
Poo, of 10°; however, the Poo, values differ 
significantly between the three watersheds (see 
summary values in Table 1) and reflect the dif- 
ferences in overlying vegetation and in the soil 
profile development with respect to amounts of 
carbonate minerals. The Huron watershed has 
Poo, Values between 10°* and 10*°, the Che- 
boygan has P.,, values ranging from 10 to 
10°'8, and the Tahquamenon watershed has the 
widest range of P. values, from 10° to 1071. 
The Huron watershed has the longest growing 
season and more cropland than the other two 
watersheds, which can help to explain the higher 
Poo, Values. 

The variability in the (IAP/K). ie Values for 
the groundwaters from the three study water- 
sheds is shown in Figure 8. The groundwaters 
within the three watersheds are at or near satura- 
tion for both calcite and dolomite. As shown in 
Table 1, calcite is on average ~1.5 times satu- 
rated in the groundwater, while dolomite is on 
average near equilibrium. Due to the solubility 
relations and the common ion effect (as shown 
in Fig. 2B), temperate soil water and groundwa- 
ter that contact both calcite and dolomite are in 
theory capable of dissolving calcite to saturation 
and then continuing to dissolve dolomite, which 
remains undersaturated. If dolomite saturation 
is attained, then calcite must become supersatu- 
rated in the solution. 

Because the groundwaters in our study are 
~1.5 times supersaturated with respect to cal- 
cite and slightly undersaturated with respect to 
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dolomite, this scenario requires that the disso- 
lution kinetics of dolomite in groundwaters be 
more rapid than the precipitation kinetics of cal- 
cite. Near equilibrium for carbonate minerals, 
reaction rates of dissolution versus precipitation 
have very different kinetic rate laws such that 
dissolution is much more rapid than is precipita- 
tion at an equivalent departure from equilibrium 
(e.g., Walter, 1985; Blum et al., 1990). Thus, 
using saturation state averages for groundwaters 
in Table 1, it is not unreasonable to expect that 


0.32 0.6 
40 
[| Huron 
a0 | Cheboygan 
# of 
samples 


Tahquamenon 


dolomite dissolution rates at a 42 IAP/K,,, 10 
value of 0.8 would be more rapid than calcite 
back precipitation rates in the same solution at a 
TAP/K. ice Value of 1.3. 

The saturation states with respect to calcite 
and calculated P,,, values for both groundwaters 
and related surface waters for the Tahquamenon, 
Cheboygan, and Huron watersheds are shown 
in Figures 9A, 9B, and 9C, respectively. The 
Tahquamenon surface waters have a wide range 
of P.o, values, from 10°! to less than 10°°. 
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Figure 8. Ranges of groundwater calcite saturation states for each of the study three water- 
sheds, shown as histograms. All three watersheds tend to have groundwaters that are slightly 
to significantly supersaturated for calcite and, thus, approach equilibrium with respect to 


more soluble dolomite. 
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These extremely low values are equivalent to 
the P.,, of the atmosphere. The average ground- 
water Pag, values differ significantly among the 
three watersheds (Tahquamenon: 10° *°%%; 
Cheboygan: 107*°*°%°; Huron: 10713#°!), 
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Figure 9. Plot of IAP/K,,,,, versus Pro, for 
surface waters relative to groundwaters for 
the three study watersheds. The Tahquame- 
non watershed (A) exhibits little difference in 
saturation state and P.o, for the groundwa- 
ters and surface waters. The Cheboygan sur- 
face waters (B) indicate a departure from the 
groundwater values and are approximately 
eight times saturated with respect to calcite. 
The Huron surface waters (C) are nearly 30 
times saturated with respect to calcite. 
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The groundwaters are saturated with calcite at 
Poo, Values much higher than that of the atmo- 
sphere. As this groundwater discharges to the sur- 
face, the gradient between the groundwater and 
the atmosphere facilitates the degassing of CO, 
from the water. Stream residence time and bio- 
logic activity also influence the loss of CO, from 
the stream water. This creates the conditions nec- 
essary for the increased saturation of calcite in 
the surface waters. The Huron watershed has the 
sharpest gradient between the Pro, value of the 
atmosphere and the groundwaters (see Fig. 9C). 
It also has the greatest degree of supersaturation 
for calcite in the surface waters. The Cheboy- 
gan watershed also has surface waters that are 
supersaturated for calcite, although not as highly 
supersaturated as the Huron (see Fig. 9B). The 
Tahquamenon watershed has only a few surface 
waters that are more saturated than the ground- 
waters; these are the waters that have lost Co, 
relative to the groundwater (see Fig. 9A). The 
degree of supersaturation with respect to calcite 
suggests that there may be losses of Ca’* and 
HCO, to carbonate reprecipitation along the 
stream drainage system. 


Mass Balancing Dissolution and Precipitation: 
Mg”/HCO, as Conservative Tracer 


In principle, it should be possible to utilize 
the concentration relations between Ca?* and 
HCO, to examine losses due to calcite precipi- 
tation within stream reaches. However, track- 
ing these changes in stream water chemistry is 
hampered by anthropogenic inputs to the sur- 
face water system. As discussed earlier, pollu- 
tion from road salts adds excess Ca**, Na*, and 
CI, but not Mg** to surface waters, and this is 
particularly significant for the Huron watershed. 
The excess Ca** makes it an unreliable tracer of 
carbonate precipitation/dissolution. Instead, it 
is possible to utilize the amounts of Mg?* and 
HCO,, which are not a significant component of 
the anthropogenic salt load, to mass balance the 
amount of dissolution versus reprecipitation of 
carbonates in the surface water system. Impor- 
tantly, only low-Mg calcite can precipitate if the 
saturation state increases via CO, loss, because 
dolomite precipitation is kinetically inhibited 
(e.g., Suarez, 1983; Kempe and Emeis, 1985; 
Komor, 1994). Thus, the Mg*/HCO, ratio of 
surface waters relative to groundwaters provides 
a means of determining whether HCO, has 
been lost via CaCO, precipitation (Mg*/HCO, 
increases) or whether the water has simply been 
diluted by rainfall or concentrated by evapora- 
tion (Mg*/HCO,- remains constant). 

The relationships between concentrations of 
Mg* and HCO, in groundwaters and surface 
waters of the three watersheds are shown in 


Figure 10. The groundwaters have a fairly con- 
stant Mg**/HCO, average ratio of 0.18, with a 
standard deviation of +0.03 (see also Table 1 
summary). This ratio proves to be a very useful 
quantity because it is determined in the ground- 
water, which is the major contributor of these 
ions to the streams. Because relatively small 
amounts of Mg” are removed during the precip- 
itation of low-Mg calcite, the Mg*/HCO, ratio 
can be used to compare differences between the 
initial groundwater chemistry and that of the 
more highly supersaturated surface waters. 

The Tahquamenon and the Cheboygan sur- 
face waters all have the same Mg*/HCO, molar 
ratio as the local shallow groundwater (see 
Fig. 10A). The relative constancy of the Mg”*/ 
HCO, ratios of groundwaters and streams in 
these watersheds suggests that processes occur- 
ring within the streams do not significantly alter 
the early imprint of weathering zone values 
on carbonate weathering fluxes from dolomite 
and calcite. The surface stream waters from the 
Huron watershed have Mg*/HCO, molar ratios 
(0.23 + 0.05) that are higher than the groundwa- 
ter ratio (see Tables 1 and 2; Fig. 10B). The rela- 
tive loss of HCO, can be attributed to localized 
precipitation of calcium carbonate within lake 
and wetland systems within the Huron water- 
shed. A more detailed discussion of the sites of 
carbonate precipitation in the Huron from the 
headwater to downstream reaches is provided in 
Szramek and Walter (2004). Importantly, the net 
effect of calcite precipitation on riverine carbon- 
ate weathering fluxes from the Huron is small 
because precipitation is only evident at lowest 
discharge conditions (Szramek et al., 2007). 


Comparisons with North American 
Drainage: Divalent Cation and Inorganic 
Carbon Budgets 


The results obtained from our study of surface 
water and groundwater geochemistry in three 
diverse watersheds established on glaciated 
landscapes of the mid-continent region can be 
placed in a larger context through comparison 
to their larger drainage systems. The two most 
significant rivers draining the North American 
continent are the Mississippi and St. Lawrence 
Rivers. As shown on the drainage basin map for 
these rivers (see Fig. 1), a significant amount 
of the drained region has been recently glaci- 
ated and is largely in contact with glacial drift 
deposits. Average composition and annual water 
and elemental fluxes were obtained for the Mis- 
sissipp1 and St. Lawrence River basins. Data 
for each of the main reaches of the Mississippi 
(Missouri, Upper Mississippi, Ohio, Tennessee, 
Arkansas, and Red Rivers) and for the three 
Michigan study watersheds (Tahquamenon, 


Geological Society of America Bulletin, May/June 2007 


Carbonate geochemistry in temperate glaciated watersheds 


HCO3 (meq/L) 


y=4.96x+ 0.25 
P= 0.84 


@ Tahquamenon 
O Cheboygan 
@ Huron 


HCOx3 (meq/L) 


Groundwater 


0 0.2 


0.4 0.6 0.8 


Mg** (mM) 


1.0 


Figure 10. Mg” versus HCO, concentrations for groundwaters (gw) (A) and surface stream (B) waters in the three study catchments. 
Groundwaters from all three watersheds have a similar Mg”*/HCO, ratio value, as shown by the regression equation. Both the Cheboygan 
and Tahquamenon watersheds have Mg”*/HCO, ratios in surface waters comparable to their respective shallow groundwaters. The Huron 
watershed has a ratio that has increased relative to the source groundwater, indicating likely back precipitation of low-Mg calcite in the 


lower reaches of the Huron River. 


Cheboygan, and Huron) were obtained from 
the USGS-National Stream Quality Accounting 
Network (NASQAN) database (Alexander et al., 
1998) using compiled data on average river dis- 
charge (in L yr'), water chemistry (Mg**, Ca’*, 
HCO, concentrations), and drainage basin area 
(km?). Data on average discharge and chemical 
composition for the mouth of the St. Lawrence 
River were obtained from Hélie et al. (2002) and 
Rondeau et al. (2005). 

The relations between Ca** and Mg** concen- 
trations for the various river drainages are shown 
in Figure 11 along with the value for world aver- 
age river water (data from Berner and Berner, 
1996). All of these North American river drain- 
ages have higher Mg** and Ca** concentrations 
than world average river water. The Mg”*/Ca** 
ratios in all these rivers fall within a range of val- 
ues between 0.30 and 0.75, indicating that dolo- 
mite and calcite weathering contribute subequal 
proportions of these divalent cations. The bud- 
gets for water, HCO,, and Mg” annual fluxes 
from the main North American river drainages 
are shown in Figure 12. The contributions of 
watersheds that fall within the recently glaciated 
region of North America (St. Lawrence, Upper 
Mississippi, Ohio, and Missouri Rivers) show 
that these are disproportionately large for both 


silicate mineral weathering rates in geochemi- 
cal models (e.g., Berner et al., 1983; Bluth and 
Kump, 1994). The comparisons made here show 
that contributions from carbonate mineral weath- 
ering, and especially the relative proportions of 
dolomite relative to calcite weathering, can actu- 
ally increase with decreasing temperature. 


SUMMARY AND CONCLUSIONS 


The suite of geochemical analyses obtained 
on the surface waters and shallow groundwa- 
ters of the Tahquamenon, the Cheboygan, and 
the Huron watersheds in the U.S. mid-continent 
demonstrates a remarkable regularity in the 
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Mg” and HCO, budgets. 

Taken together, these comparisons show that 
dolomite and calcite weathering are both impor- 
tant contributors to riverine solute budgets and 
that temperate, recently glaciated watersheds 
are disproportionately large contributors to car- 
bonate weathering fluxes. Increased global tem- 
perature is widely considered to be a control on 


Geological Society of America Bulletin, May/June 2007 


0.50 0.75 


1.0 1.25 


Ca** (mM) 


Figure 11. Average Mg”* versus Ca” concentrations for the three study watersheds in Michi- 


gan (Tahquamenon, Cheboygan, and Huron) 


versus average values for the mouth of the 


St. Lawrence River and Mississippi River tributaries (Ohio, Upper Mississippi, Missouri, 
Ohio, Tennessee, Red, and Arkansas Rivers). The composition of world average river water 


is also indicated. 
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N. American Drainage Budgets 
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Figure 12. Annual drainage budgets for the St. Lawrence and tributaries of the Mississippi 
River for water, HCO,, and Mg” fluxes. Tributaries within the extent of Pleistocene glacia- 
tion (see Fig. 1) are indicated by shaded bars in each histogram. The numbers within the 
histograms for HCO, and Mg” represent the average concentration of the respective species 
in each of the tributaries. The trends indicate that most of the Mg** and HCO, fluxes from the 
main North American drainage systems are derived from temperate, glaciated watersheds. 


parameters controlling the dissolved inorganic 
carbon mass balance within these watersheds. 
Despite their large spatial separation and dif- 
ferent glacial histories and types of glacial 
drift, groundwaters from all three watersheds 
in Michigan have a very consistent Mg”*/ 
HCO, mole ratio value of ~0.2. The differences 
in average HCO, value of shallow groundwa- 
ters are controlled by equilibrium with respect 
to different P., values. Surface water Mg** 
and HCO, concentrations for the Cheboygan 
watershed are very close to those of the shallow 
groundwaters, while those from Tahquamenon 
and Huron watersheds tend to be lower than 
their respective groundwaters, likely because of 
groundwater dilution by runoff. 

Because of the local use of CaCl, de-icing 
salts (Jones and Sroka, 1997), Ca** concentra- 
tions are poor indicators for carbonate dissolu- 
tion or precipitation mass balances. Instead, the 
relatively constant value of the Mg*/HCO, 
ratio of groundwaters can be used to evaluate 
changes in surface waters due to calcite back 
precipitation. The only surface streams that 
have markedly different Mg*/HCO, ratios 
from groundwaters are in the Huron watershed, 
where HCO, appears to have been selectively 
removed by carbonate precipitation. However, 
even in the Huron watershed, the net effect 
of calcite precipitation on riverine carbonate 
weathering fluxes is small because calcite pre- 
cipitation is only evident at lowest discharge 
conditions (Szramek et al., 2007). 

The degree of calcite saturation and P,. val- 
ues calculated for surface waters and groundwa- 
ters also show significant differences among the 
three watersheds. Both the Cheboygan and the 
Huron watersheds have groundwaters that are 
~1.6 times saturated with respect to calcite and 
near equilibrium with respect to dolomite, at dif- 
ferent subsurface P.,. values. Stream and river 
waters in both watersheds have evolved to higher 
calcite saturation states. However, the Huron 
River has evolved to a much higher degree of 
calcite supersaturation, generally between 8 and 
10 times supersaturated. The range of Mg**/Ca** 
values observed in groundwaters and surface 
waters of the three study watersheds in Michi- 
gan suggests that dolomite and calcite dissolu- 
tion contribute subequally. 

We suggest that the increase in the relative 
and absolute solubility of dolomite relative 
to calcite at temperatures below 30 °C plays 
an important role in regulating the weather- 
ing fluxes of these two minerals. Comparisons 
among the main North American river drainages 
support this contention, showing that both Mg** 
and HCO, riverine concentrations and annual 
fluxes are heavily weighted toward drainages 
located in the temperate, glaciated portion of the 
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continent. The lack of significant modification 
of inorganic carbon fluxes by back precipitation 
in surface waters is also an important character- 
istic of these watersheds. 

These results have significance in terms of 
the short- and long-term carbon cycle and the 
role played by carbonate mineral weathering. In 
the short term, the likely increases in the size 
of the terrestrial and atmospheric carbon pools 
due to fossil-fuel burning would be expected to 
increase carbon cycling rates in soils and, hence, 
weathering zone P,, 0, values. Carbonate mineral 
weathering fluxes to surface waters and shallow 
groundwaters should increase, with influence 
on both surface water chemistries and on global 
fluxes to the ocean reservoir. 

Although carbonate mineral weathering is 
generally considered unimportant in long-term 
carbon cycle models, the significant role of 
dolomite weathering suggested here may require 
revision of that assumption. Dolomite is widely 
known to have been much more common as a 
primary marine mineral in other geologic time 
periods (e.g., Holland and Zimmermann, 2000). 
Today, dolomite weathering on continents pro- 
duces HCO, that cannot be balanced by dolo- 
mite precipitation in modern oceans. Thus, 
asymmetry of the Ca** and Mg** sources during 
weathering on continents (dolomite and calcite 
weathering) and the fate of the resultant HCO, 
in the marine environment (calcite versus dolo- 
mite precipitation) merit additional study. 
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